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We study a model of n interacting fermions in a disordered potential, which is assumed to 
generate uniformly fluctuating interaction matrix elements. We show that the ground state magne- 
tization is systematically decreased by off-diagonal fluctuations of the interaction matrix elements. 
This effect is neglected in the Stoner picture of itinerant ferromagnetism in which the ground-state 
magnetization is simply determined by the balance between ferromagnetic exchange and kinetic 
energy, and increasing the interaction strength always favors ferromagnetism. The physical origin of 
the demagnetizing effect of interaction fluctuations is the larger number K of final states available 
for interaction-induced scattering in the lower spin sectors of the Hilbert space. We analyze the 
energetic role played by these fluctuations in the limits of small and large interaction U . In the 
small U limit we do second-order perturbation theory and identify explicitly transitions which are 
allowed for minimal spin and forbidden for higher spin. These transitions then on average lower 
the energy of the minimal spin ground state with respect to higher spin; we analytically evaluate 
the size of this reduction and find it to give a contribution A s oc nU 2 /A to the spin gap between 
the two lowest spin ground-states. In term of an average effective Hamiltonian, these contributions 
induce a nU 2 S 2 /A term which decreases the strength of the ferromagnetic exchange, thereby delay- 
ing the onset of Stoner ferromagnetism, and generate a second, larger S term oc S 3 which results in 
a saturation of the ground-state spin before full polarization is achieved, in contrast to the Stoner 
scenario. For large interactions U we amplify on our earlier work [Phys. Rev. Lett. 84, 3938 (2000)] 
which showed that the broadening of the many-body density of states is proportional to \fKU~ and 
hence favors minimal spin. Numerical results are presented in both limits. After evaluating the 
effect of fluctuations, we discuss the competition between fluctuations plus kinetic energy and the 
exchange energy. We finally present numerical results for specific microscopic models and relate 
them to our generic model of fluctuations. We discuss the different physical situations to which such 
models may correspond, the importance of interaction fluctuations, and hence the relevance of ous 
results to these situations and recall an experimental setup which we proposed in an earlier worku 
to measure the importance of interaction fluctuations on the ground-state spin of lateral quantum 
dots in the Coulomb blockade regime. 

PACS numbers : 73.23.-b, 71.10.-w, 75.10.Lp 

I. INTRODUCTION 

A. Stoner effect and disorder 

More than fifty years ago Stoner proposed a simple route to,ferromagnetism in itinerant systems based on the 
competition between one-body and interaction (exchange) energyH. The repulsive interaction energy can be minimized 
when the fermionic antisymmetry requirement is satisfied by the spatial wavefunction, as the overlap between different 
wavefunctions is then minimal. This effect favors the alignment of spins and, if the interaction is sufficiently strong, 
results in a large ground-state spin magnetization. This mechanism is the primary origin of Hund's first rule in atomic 
physics. In contrast, when the interaction is weak minimal spin is favored, since in order to align spins electrons must 
be promoted from lower, doubly-occupied levels to higher singly-occupied levels and the cost in one-body energy is 
prohibitive. Because the Pwt principle is essentially local, ferromagnetism in metals has been studied within models 
such as the Hubbard modeltra which only retain the short-range part of the electronic interaction; the long-range part 
of the interaction being assumed to give spin-independent contributions to the ground-state energy (the capacitance 
or charging energy). In the case of a Hubbard interaction, only pairs of electrons of opposite spin interact. The 
number of such pairs is a monotonically decreasing function of the total magnetization ~ [(n/2) 2 — a 2 ] where n is the 
number of electrons and a the total spind. On the other hand, as just noted, flipping a spin requires the promotion 
of an electron to a higher one-body level and in the case of a finite system with a discrete spectrum of average 
spacing A, a magnetization a requires an energy c 2 A. A simple first order perturbation treatment shows then that a 
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sufficiently strong interaction results in a finite magnetization, when the corresponding reduction in interaction energy 
counterbalances the increase in kinetic (one-body) energy 

(A-V c )o- 2 = (1.1) 

This occurs when the typical exchange interaction V c between two states close to the Fermi energy is equal to the 
one-particle level spacing which for a Hubbard interaction U(f,f l ) = US(r — f) reads 



V C = U C J dr\ipa(r)\ 2 \'ipp(r)\ 2 = A (1.2) 

The upper bar indicates an average over wavefunctions in the vicinity of the Fermi level. In a clean system this gives 
U c = A and this threshold is known as the Stoner instability. As both the kinetic energy and the interaction energy 
have the same parametric dependence on the magnetization a, reaching this threshold results in a second order phase 
transition to a ferromagnetic phase, the divergence of the magnetic susceptibility, and a macroscopic magnetization. 

Quite naturally one may wonder in what way does the presence of a disordered potential modify this Stoner picture, 
and this question has recently attracted a lot of attention, both in the context of bulk metals (i.e. infinitely extended 
systems with diffusive eigenstates) and in finite-sized metallic systems such as quantum dots and metallic nanoparticles. 
Two types of questions have been considered: 1) The effect of a disordered potential on the average threshold for 
the Stoner instability. 2) The statistical properties of the threshold in an ensemble of mesoscopic metallic samples. 
Both aspects have been recently investigated theoretically. For the bulk case, it has been known for some timeQ that 
within perturbation theory disorder enhances the exchange effect in the susceptibility; recently Andreev and Kamenev 
constructed a mean-field theory which they argue describes the Stoner transition]] and found a significant reduction 
of the Stoner threshold in low-dimensional disordered systems due to correlations in diffusive wavefunctions which 
enhance the average exchange term. In the framework of the same mean-field approach which neglects the fluctuations 
of the interactions, but takes into account those of the one-body spectrum, Kurland, Aleiner and Altshulcr proposed 
that below, but in the immediate vicinity of the Stoner instability, there is a broad distribution-of magnetization and 
that eachisample's free energy is characterized by a large number of local magnetization minimac3. Brouwer, Oreg and 
Halperintl considered the effect of mesoscopic wavefunction fluctuations on the exchange interaction and found that 
their effect was to increase substantially the probability of nap-zero spin magnetization in the ground state before the 
Stoner threshold is reached. Baranger, Ullmo and GlazmanB suggested that the observed "kinks" in the parametric 
variations of Coulomb blockade peak positions (e.g. as one varies an external magnetic field) could reflect changes 
in the ground state spin of the quantum dot. It was noted that the statistical occurence of nonzero ground-state 
magnetizations can account for the absence of bimodality of the conductance peak spacings distribution for tunneling 
experiments with quantum dots in the Coulomb blockade regimeE3~li3. Another aspect of large disordered metallic 
samples is that the Stoner threshold can be locally exceeded, while the exchange averaged over the full sample has a 
value well below the threshold. In this case one may expect that localized regions with nonzero magnetization will 
be formed even though the full system is nonmagnetic. This scenario has been investigated by Narozhny, Aleiner and 
LarkinEj who also considered the effect of such local spin droplets on dephasing. They found that the probability to 
form a local spin droplet though exponentially small, does not rigorously vanish as it would in a clean system, and that 
neither this probability, nor the corresponding spin depends on the droplet's size. In a different approach focusing 
on the large interaction regime close to half-filling, Eisenberg and Berkovits numerically found that the presence 
of disorder may stabilize Nagaoka-like ferromagnetic phases at larger number of holes (> 2)E3. Finally, Stopa has 
suggested that scarring of one-body wavefunctions in a chaotic confining potential may lead to strong enhancements of 
the exchange interaction and to the occurence of few-electron polarization in finite-sized systemstj. Thus the general 
message of these works is that disorder tends to favor novel magnetic states over paramagnetic states. 



B. Overview and outline 

In a recent letterB, we pointed out a competing effect of interactions in disordered systems which reduces the 
probability of ground-state magnetization and hence favors paramagnetism. This effect had not (to our knowledge) 
been treated in any of the previous works on itinerant magnetization of disordered systems. The wojjka cited above 
neglect the effect of disorder in inducing fluctuations in the off-diagonal inter actio n matrix elementsBEj. However it 
is well-known from studies of complex few-body systems like nuclei and atomsLZMia that the band-width of the many- 
body density of states in finite interacting fermi systems is strongly modified by the fluctuations of these off-diagonal 
matrix elements already at moderate strength of the interactions. Such studies did not directly address the effect of 
this broadening on the ground-state spin of the system. However our extension of these models immediately revealed 
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that these fluctuations are largest for the states of minimal spin, due to the larger number of final states (non-zero 
interaction matrix elements) for interaction-induced transitions (we will review this argument below) . This effect then 
significantly increases the probability that the extremaLilow-energy) states in the band are those of minimal spin 
and opposes the exchange effect. In our earlier workso'Ej we focused on the regime of large fluctuations to deduce 
the scaling properties of the ground-state energy as a function of spin, and verified these scaling laws with numerical 
tests. In the present work we will review and extend these results for large fluctuations, but we will focus mostly on 
the perturbative regime of small U. While in this regime the correction to the ground-state energy due to fluctuations 
is small by assumption, one is able to evaluate these corrections analytically and show that they favor minimal 
spin for an arbitrary number of particles. Specifically, the larger number of interaction-induced transitions for lower 
spin leads to more and larger interaction contributions to the (negative) second-order correctkinJo the ground-state 
energy in each spin block. This illustrates explicitly the "phase-space" argument introduced mllrEJj which implies that 
fluctuations generically suppress magnetization. We expect this effect to be significant in quantum dots where it will 
reduce the probability of high spin ground states. We recall that the ground state spin of lateral quantum dots can be 
experimentally-determined by following the motion of Coulomb blockade conductance peaks as an in-plane magnetic 
field is applied^. Therefore the strength of the demagnetizing effect of fluctuations of interaction is experimentally 
accessible. 

The paper is organized as follows. In section || we start by an explicit derivation of our model and describe its main 
features. In section III we begin for pedagogical reasons with an analytical treatment of the model for the case of only 
two particles, in both the perturbative regime of weak off-diagonal fluctuations and the asymptotic regime where they 
dominate. Section IV will be devoted to a second order perturbative treatment of the model for an arbitrary number 
of particles; this will be followed in section [y] by a discussion of the magnetization properties of the system's groi 
state in the asymptotic regime. As noted above, some of the results presented there have already been presented ir 
but are nevertheless included to make the article self-contained. In the next section VI we consider the competition 
between exchange and fluctuations in more details, both from the point of view of average Stoner threshold and in 
terms of probability of finding a polarized ground-state. We will see in particular that the off-diagonal fluctuations 
induce a term ~ a 2 in the Hamiltonian which delays the Stoner instability and a second term ~ a 3 which strongly 



suppresses the occurence of large ground-state spins even above the Stoner instability. In section VII we consider 
more standard microscopic models for disordered interacting fermions and relate their properties to our generic model 
of fluctuating interactions. We determine the conditions to be satisfied in order for the results obtained from our 
random interaction model to be relevant in different physical situations. Finally we summarize our findings, put them 
in perspective and discuss possible extensions of this work in the final section VIII . 



II. DERIVATION OF THE MODEL 



Our starting point is a lattice model for fermions in a disordered potential coupled by a two-body, spin-independent 
interaction of arbitrary range. We make a unitary transformation to the basis of single-particle eigenstates of the 
disordered potential and introduce the assumption that the single-particle states are random and uncorrelated. Upon 
averaging over disorder we arrive at a completely generic model describing both the nonvanishing average interactions 
(exchange, charging and BCS) and the statistical fluctuations in both the one and two-body terms. Finally we introduce 
the assumption that all interaction matrix elements have the same statistical variance. Hence our construction excludes 
both one-body integrable systems and strongly-localized systems. We also note that with this assumption geometric or 
commensurability effects (such as spin waves or antiferromagnetic instabilities) cannot be captured by our treatment, 
as the statistical character of the construction erases most real-space details of the model. 

We consider the following tight-binding Hamiltonian for n spin-1/2 particles 

n = n + u = Y^ K j 4 A* + U ^ ~ fi n%n i + H u(Q)n iA n iA 

i,j;s i^ij i 

s =t,l is a spin index, d\ (di,a) creates (destroys) a fermion on the i th site of a D-dimensional lattice of linear 
dimension La and volume fi = (La) D . This latter quantity defines the number m/2 = fl/a D of spin-degenerate one- 
body eigenenergies which we will refer to as orbitals in what follows, a = 1 is the lattice constant. H.q is a one-body, 
spin independent, disordered Hamiltonian with eigenvalues e a and eigenvectors ip a , i.e. one has 

Hol^a) = eal^a) = e a y]i(i a (i)\i) (2.1) 

i 

where \i) refers to a lattice site ket. 
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We assume that Tio has no degeneracy besides twofold spin-degeneracy and distribute the m/2 different one-body 
energies as e Q € [0;m/2] so as to fix A = 1 without spin degeneracy. Below we will discuss three different eigenvalue 
distributions : constant spacing distribution^!] (e Q = (a — 1)A, note that due to the level degeneracy the single- 
particle level spacing is A, whereas A/2 is the mean level spacing), randomly distributed e a with a Poisson spacing 
distribution, or with a Wigner-Dyson spacing distribution. Finally I4(i—j) is the electron-electron interaction potential 
and rii = ^ s nj lS = ^2 S d\ s di tS . The Hamiltonian is Spin Rotational Symmetric (SRS) so that both the total spin \S\ 
and its projection S z commute with the Hamiltonian and the corresponding eigenvalues a and o~ z are good quantum 
numbers. This results in a block structure of the Hamiltonian TL which will be described in detail below. Performing 
the unitary transformation defined by 

5>a(i)4,, = <•> (2.2) 

a 

we rewrite the Hamiltonian as 

W = £a4,s C <»i« + ^Q,/3 C a,s C ks' C|5 ' s ' C ')'> s (2-3) 

where the Interaction Matrix Elements (IME) are given by 

These IME induce transitions between many-body states differing by at most two one-body occupation numbers. 
The distribution and properties of the IME depend on both the range of the interaction potential and the one-particle 
dynamics. If there are conserved quantities other than energy in the one-particle dynamics (and hence good quantum 
numbers describing the one-body states) this will lead to selection rules in the IME; the extreme case of this would 
be an integrable one-particle hamiltonian for which a complete set of quantum numbers exists. Selection rules greatly 
reduce the number of allowed interaction-induced transitions, and lead to a very singular distribution of IME (this 
is most easily seen by considering a clean hypercubic lattice model with Hubbard interaction). Perturbing a clean 
lattice with a disordered potential destroys translational symmetry and these selection rules disappear, which induces 
a crossover of the distribution of IME from a set of S— functions to a smooth distribution. In Figffl we illustrate this by 
plotting the distribution of IME for a one-dimensional lattice model with pa-site disorder, nearest and next-nearest- 
neighbor hopping and a Hubbard interaction as described e.g. in referenceEil. 

The key assumption of our model is that such a smooth distribution of interaction matrix elements exists and that 
in fact all matrix elements which preserve SRS have the same nonzero variance (these matrix elements may vanish on 
average of course). This assumption rules out both the case of integrable one-body dynamics as discussed above, and 
the case of strongly-localized wavefunctions for which interaction matrix elements between states separated spatially 
by more than a localization length will have different (and much smaller) variance than those in the same localization 
volume. Our assumption is reasonable for metallic disordered states with a randomness generated by either impurities 
or chaotic boundary scattering. 

With this motivation, we assume that the fluctuations of the off-diagonal U^'i are random with a zero-centered 
gaussian distribution of width U. Only matrix elements U^'g, U^'g and U@>P have nonzero averages {U"'g), {U^'g) and 
(Ua'a) which lead (respectively) to mean-field charge-charge, spin-spin and BCS-like interaction terms. Note that the 
average of both the exchange and BCS terms is dominated by the short-range part of the interaction, and that IAbcs 
vanishes as time- reversal symmetry is broken. Consequently, the electronic interactions give us four contributions. 
The first three are the average charge-charge, ferromagnetic spin-spin and BCS terms that we just discussed and 
which can be written as (n a = s c a,s c a,s and n — "<>•) 



lA aV q = Ucc + Uss + l^BCS = 



n(n + l)/2 -XUS-S + Y,< U ™) c *A,l c P.l c M ^ 



where we have introduced spin operators S a = (1/2) J2 S t c a s ^s,tc a ,t and S — J2 a Note that the strength of the 
average ferromagnetic exchange term has been written in units of the rms fluctuation U, i.e. we have introduced a 
parameter A which is the ratio of the average exchange to the fluctuations 



X = 2(U^)/U (2.6) 
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much in the same spirit as the usual Stoner picture where another energy ratio (U^'p) / A, between the exchange 
energy and the one-body energy spacing at the Fermi level, is the relevant parameter. 

The fourth interaction contribution to our model hamiltonian goes beyond the mean-field approximation and 
contains the off-diagonal fluctuations of the electronic interactions 

u f= E J2^iA, S 'C S , s ^, s (2.7) 

a,/3;7,<5 s.s' 

Having removed the average interactions, we now assume that both the diagonal and off-diagonal IME U^'p have 

zero-centered uncorrelated gaussian distributions PiU^'s) oc e 1 a <P' 1 of width U. We stress that in general, not 
all IME have the same variance but being interested in generic features of the interaction, we will neglect these variance 
discrepancies. Uf contains three kind of matrix elements, the variances of which depend on the number of transferred 
one-body occupancies between the connected Slater determinants. Diagonal matrix elements Uf' 1 — (I\Uf\I) 
denotes a Slater determinant) have a variance ~ n(n — l)U 2 /2, one-body off-diagonal elements that change only one 
occupancy have a variance ~ (n — 1)U 2 whereas generic two-body off-diagonal matrix elements inducing transitions 
between Slater determinants differing by exactly two occupancies, have the generic variance U 2 . In diagrammatic 
language, these discrepancies occur due to the presence of up to two closed loops in the diagram corresponding to 
these matrix elements, each loop corresponding to a sum over 0(n) uncorrelated IME. 
Our full Hamiltonian then reads 



H = H +U avg +Uf (2.8) 
The mean-field Hamiltonian proposed inEl was constructed along similar lines but neglects the fluctuations of 



interaction Uf and is thus embedded in the above Hamiltonian (2.8). Consequently, all results derived there can be 
obtained from the treatment to be presented below after setting the strength of fluctuations U ^ {). In a condensed 
matter context this is justified in the limit of large .conductance g — > oo. As recent experiments in quantum dots 
seem to be consistant with a conductance g w 6 — 8EJ, it is a priori not obvious that Uf can be neglected. We also 
stress that both the Random Matrix Theory (RMT) symmetry under orthogonal (or unitary) basis transformation in 
the one-body Hilbert space (which in metallic samples is satisfied for energy scales smaller than the Thouless energy 
E c = gA\H), and the SU(2) symmetry under rotation in spin space are satisfied by each of the three terms in the 
above Hamiltonian. 

The charge-charge mean-field contribution results in a constant energy shift of the full spectrum and has thus no 
influence on the ground-state spin; we therefore neglect it henceforth. This must however be kept in mind, as it is 
for instance well known that including self-consistently the mean-field charge-charge contribution of the interactions 
(e.g. | -.in |[ a 1 Hartree-Fock approach) leads to significant corrections to the one-particle density of states at the Fermi 
leveEaEl The BCS term gives rise to superconducting fluctuations for a negative effective interaction in the Cooper 
channel (U^'P) < 0. We shall only consider disordered metallic samples which have (C/^'f ) > 0. In this case the 
renormalization group flow brings the BCS coupling to zerocil. We thus also neglect this term and set U avg = IA S8 . 
Note however that the presence of a nonzero (repulsive or attractive) BCS coupling may stabilize a paramagnetic 
phase. 

After these considerations we reach our model Hamiltonian 

ck a,/3;7,<5 s,s' 



Due to the SRS that we imposed on the original Hamiltonian ( |2.l| ), the interaction commutes with the total 
magnetization | *S" | 2 and its projection S z so that the Hamiltonian acquires a block structure where blocks are labelled 
by a quantum number a z and subblocks of given a > \a z \ appear within each of these blocks. Each block's size is 

given in term of binomial coefficients as N(a z ) — (^^-a- ) (n/2+<r )' while the size of a subblock of given a is given 
by N(<r) = N(a z = a) — N(a z = a + 1). Due to SRS it is sufficient to study the block with lowest projection a z = 
(1/2) for even (odd) number of particles, as all values of a will be included in this block. For simplicity, we will 
consider an even number n of particles in the initial discussion presented below, and will generalize the discussion 
later on to include odd n, highlighting the main differences between the two cases. It is important to remark that 
both a and a z are not only good quantum numbers for the full Hamiltonian, but also individually for Ti.Q, lA avg and 
Uf. This allows us to consider each of these terms separately and in the next two chapters we will make use of this 
property, first neglecting U ss : as it only generates constant energy shifts within each sector, it can be added after the 
restricted problem Ho + Uf has been solved. 
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In ( |2.7| ) the sums in both the spin and orbital indices are not restricted, i.e s,s' =1,1 and a, 0,^,5 = 1,2, ...m/2. 
It is both convenient and instructive to rewrite it as 

u f= E E Kpl {s z )T 1>s { Sz ) + l - y, (2-io) 

a>,3;7<<5 s z =0,±l a>,3:7<<5 

where we have introduced totally symmetric and antisymmetric matrix elements 

v2i = uli + ufc - u s a } - uj;l (2.11) 

as well as two-body creation and destruction operators for either singlet-paired 

S a ,/3 = ( C 1,T C /3,1 _ C a,l C /3,t)/^ > ^a,a = c l,T C Ll ( 2 -12) 

or triplet-paired fermions 

^(0) = (4, T 4 a +4 a ct )T )/V2 ; T^(s) = cl s cl a , S =U (2.13) 

As we consider fully uncorrelated IME U^'g, both the symmetrized and antisymmetrized matrix elements have the 
same variance which for no doubly appearing indices reads 

- 2 (W$) = ° 2 (K#) = 4C/ 2 (2.14) 

In principle, the ratio of the variances strongly depends on microscopic details, in particular the range of the interac- 
tion. For instance, it can easily be seen that a 1 (v^'^j I® 2 y^a p) e lP> 1] ana - that ^ ac ra ti° vanishes for a Hubbard 
interaction. We will neglect this discrepancy however, but note that anjncreased variance of the symmetrized IME 
with respect to the antisymmetrized ones favors a low spin ground-stateEi 

The hamiltonian can now be regarded as acting on singlet or triplet bonds between levels. SRS is then reflected 
in the simple statement that the destruction of a bond between two fermions must be followed by the re-creation of 



a bond of the same nature. We note that the triplet operators ( 2.13 ) create either a a z = 0, a = 1 or a a z = ±1, 
a = 1 two-fermion state in a fixed spin basis. A rota tion in spin space would bring the operators in ( 2.13| ) into one 



another and the first three terms in the brackets in (2.10) are not individually SRS but must be considered as one 
single spin-conserving operator. We illustrate this point in Appendix A|, where we e valu ate th e effe ct of this operator 
acting on a four-particle state with two double occupancies. Note also that from (2.4) and ( [2.11 ), a purely on-site 
interaction influences only the sing let channel as in this case the antisymmetrized IME vanish identically. 

The procedure leading to ( 2.10| ) amounts to a projection of the interaction operator onto the two irreducible 



representations of the two-fermion symmetry group. In this way the tw o-bod y singlet matrix elements are explicitly 
separated from their triplet counterparts, and the rewriting leading to ( [2.10 ) allows us to fopnulate the many-body 



problem in term of two-particles bonds of different nature in a similar way as the authors o Any even n-fermion 
state is represented as a n/2-boson state where each boson has either spin a = or 1. These bosons can be constructed 
by acting on the vacuum |0) with an S or a T operator respectively, and the spin of these composite bosonic states 
depends on the bond between the two fermions, i.e. whether the fermionic antisymmetry is supported by the spin or 
the spatial degrees of-freedom. Alternatively, this means that for a n-body state of total spin a, the number of triplet 
bonds is given by ctcj. Also double orbital occupancies result in singlet bonds, so that their number is restricted 
to [0,n/2 — a]. This construction leads however to an overcomplete basis for n > 6. We were unable to propose a 
systematic reduction to an orthonormal set of states nor are we aware of any such systematic construction in the 
literature. For the computations to be carried below it will however be sufficient to know that such a basis can in 
principle be constructed (via e.g. reduction and orthogonalization of the constructed overcomplete basis), and how to 
construct it for the special case of four particles above the filled fermi sea, as those are the only states one encounters 
when doing second order perturbation theory for the levels of lowest energy in the a = and 1 sectors. 

Equation ( |2.1C| ) helps us see the key qualitative point of our work. In second order perturbation theory Uf will 
generate transitions in each spin sub-block between the ground state and excited states differing by two occupation 
numbers (or less). Both the triplet and singlet terms will generate transitions, but there are certain types of transitions 
which can be generated by the singlet term which cannot be generated by the triplet term. For instance the triplet 
operator cannot generate transitions to final states with additional double occupancies nor is it possible to scatter 
a triplet bonded pair into a double occupancy (see Fig. 0). As the magnetization increases, the number of singlet 
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transitions decreases accordingly as the number of singlet two-particle bonds in a many-body state obviously decreases 
with its total magnetization. Eventually, when a is maximal, only triplet transitions survive and we can readily 
conclude that the number of two-body transitions is a monotonically decreasing function of the magnetization as is 
therefore the number of (energy-decreasing) second order contributions. We will see below that this condition on 
the available volume for phase-space scattering is crucial for the ground-state magnetization properties, both in the 
perturbative regime (C//A -C 1) and in the asymptotic limit of dominant fluctuations (U/A 1). It is important 
to understand that the relevant variable here is the number of transitions and not the size of the Hilbert space; the 
block size N(a) is in general (for a sufficient number of particles) a nonmonotonic function of a, as on one hand 
N(a = n/2) < N(a = 0) (or N(a — n/2) < N(<r — 1/2) for odd number of particles), whereas on the other hand, 
and in the limit a <C n/2 <C to/2, it can be shown using Stirling's formula that dN(a)/da > 0. Except for very few 
particles, N(a) has its maximum at a finite magnetization, whereas the number of transitions is always maximum for 
a = 0. 

We close this introductory chapter with a brief historical survey of random interaction models similar to ( pTg) - |2 .10 ). 
These models originated in nuclear physics and are based on similar principles as those which led Wigner to propose 
the gaussian ensembles of random matrices, with the additional requirement that they represent particles interacting 
via a fc-body interaction. Only when the rank k of the interaction is equal to the number n of particles does one 
recover the Wigner gaussian ensembles. Physically, interactions are in principle not random per se, however once one 
postulates the invariance of the one-body Hamiltonian matrix ensemble under unitary (i.e. basis) tran sformation, 
a postulate motivated e.g. by a chaotic one-body dynamics, random IME naturally appear (see (2.4)), and this 
results in a similar invariance for the many-body Hamiltonian ensemble and the associated probability distribution 
P(H) (x exp(— Tr7i 2 /2). The first proposed model with random interactions was the fermionic Two-Body Random 
Interaction Model .(TBRIM) for spinless fermions which was introduced independently by French and Wong and 
Bohigas and FloresEl This model is essentially a spinless version oiUf. While significant deviations from the usual 
gaussian ensemble of random matrices were found in the tails of the spectrum - in particular the MBDOS for n>2 
has a gaussian, not a semicircular shape - these authors found no significant differences ja the spectral properties 
at high excitation energy. (This latter finding has been however challenged very recentlyE3 and may be due to the 
smallness of the systems considered.) More recently this spinless TBRIM was extended with a one-body part and 
it has been discovered that the critical interaction [Strength U c at which WD statistics sets in is governed by the 
energy spacings A c between directly coupled siatesE3. This model and similar ones have also been studied in the 
framework of quantum chaos in atomic physicsEj, in particular the thermalization of few-body isolated 
attracted a significant attentiontlrEj, and more recently in solid-state physics to study quasinarticle lifetime 
fluctuations of Coulomb blockade conductance peak spacings and heights in quantum dotsoLj. In a solid-state context 
however the invariance of the one-body Hamiltonian under basis transformation is satisfied only in pan energy interval 
of the order of the Thouless energy E c = gA around the Fermi energy, where g is the conductances. Wavefunctions 
correlations become stronger and stronger beyond the Thouless energy where IME start to decay algebraically as 
a function of the energy. It is thus reasonable to consider our random interaction modcL-as an effective truncated 
Hamiltonian in an energy window given by the Thouless energy E £ [Ef — E c , Ep + E c ]lB, so that the number of 
particles and orbitals behave as m,n ~ g. Nuclear shell models may also be represented by randomly interacting 
models, differing from the original TBRIM in the presence of additional quantum numbers like spin, isospin, parity 
and so fortr£3. Most of those models consider the limit of dominant fluctuations U/A 3> 1 and quite unexpectedly, 
it has been found that even in this regime, random interactions may result in an orderly behavioio, in particular a 
strong statistical bias toward a low angular momentum ground-state. In particular, for the special case of an angular 
momentum restricted to j z — ±1/2, therfltobability of finding a zero angular momentum ground-state for an even 
number pinucleons reaches almost 100%Bt3. While the reasons for this behavior in the asymptotic regime are still 
not cleaio, we will see below that a strong bias toward low angular momentum ground-state results from a stronger 
broadening of the Many-Body Density of States (MBDOS) in the low spin sector, associated with a larger number 
of off-diagonal transitions. The same phenomenon with qualitatively the same origin will be shown to influence the 
ground-state magnetization in the perturbative limit. 



III. THE CASE OF n = 2 FERMIONS 



For n = 2 particles, only the sectors a = and 1 exist whose size is given by N(a) = to/2(to/2 + 1 — 2er)/2. In each 
sector, the interaction matrix Uf is a GOE matrix (the number of particles is equal to the rank of the interaction) 
and all Hamiltonian matrix elements are non-zero and have the same variance. For simplest case of two orbitals 
(to/2 = 2) one can demonstrate the magnetization reducing effect of interaction Fluctuations by an argument which 
is exact for all values of the off-diagonal fluctuations U . The two orbitals are spin-degenerate and have energies t\ = 
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and €2 = A > 0. In the absence of interaction fluctuations, the three eigenvalues in the a = sector are 0, A and 2A. 
Switching on the interaction, the determinant of the Hamiltonian matrix in the time-reversal symmetric case can be 
written 



\w# «,2 2 ) 2 + «i) 2 (2A + \w^) + {W$f(A + \wl$) 



(3.1) 



Every single term in this expression has a symmetric distribution, i.e. an equal probability of being positive or 
negative, except for a term —((W22) 2 + 2(VF L 1 2 1 ) 2 )A/4 which is always negative. It results that the determinant has a 
higher probability of being negative which in its turn means that the lowest eigenvalue (which vanishes at U/A = 0) 
is statistically more often negative than positive when U is switched on - it is more likely to be reduced than increased 
by the off-diagonal fluctuations. Simultaneously, and in absence of exchange, the energy of the only a = I level is 
given by A + 2 /2, so that the fluctuations lower or increase it with equal probability. Hence fluctuations always 
increase the average spin gap in this case. 

We next consider an arbitrary number of orbitals, m. First consider the limit of dominant fluctuations U/A 3> 1. 
H. ~Uf is then a GOE matrix and its MBDOS is well approximated by a semicircle law (E 2 < E 2 ) 



pgoe = ^eW)\! e ^- e2 (3 - 2) 

where E (a) w y/2N(a)U. This expression is not exact however as there are corrections in the tail of the distributionEj 
as one can see on Fig. |[ These corrections behave as 0{N^ 1 / & ) while the level density there is ©(TV 1 / 6 )!!!!, i.e. the 
number of levels outside the semicircle is independent of N (and hence of m) and for simplicity we will neglect these 
corrections in what follows. 

Henceforth we shall be focusing attention on the ground state in each spin sector and the gaps between these 
ground states, so it is useful to adopt the standard term in nuclear physics for the lowest levels, of a given spin or 
angular momentum, yrast levels. In the current model, in the asymptotic regime of large fluctuations U /A 3> I , (and 
neglecting the exchange interaction) we can approximate the energy of the yrast states by Eq(<j) and hence readily 
predict that the average a = yrast energy will be lower than its a = 1 counterpart by an amount 



A s = So,a=i - £o, CT =o « U(y/m(m/2 + l)/2 - y/m(m/2 - l)/2) =U+ 0( — ) (3.3) 



m 



i.e. on average there is a spin gap for U /A 3> I in the large m limit. N ext we can calculate the average energy of the 
first excited a = level £i, CT =o via integration of the average MBDOS (|3.2[) as 



Si,a=o - £a,*=a = 0(N(a)-^ 6 ) = 0( m -^ 3 ) (3.4) 

In the relevant limit m 3> I the splitting between this first excited a = level and the a = yrast level is negligible 
and both states are below the a — 1 yrast level by a gap of order U, independent of m. This calculation can be 
extended to higher a — excited states and the result suggests that on average there is a large number (^(m 1 / 3 ) of 
cr = levels which have a lower energy than the first spin excited state. Remember however that we have neglected 
corrections to the tails of the density of states, and it turns out that these corrections result in an m-independent 
number p w 3 of a = levels in the spin gap as shown by the numerical data presented in Fig. ||. In Fig. [5] we 



show a numerical check which confirm the validity of (3.3) up to prefactors which are due to additional correlations 



between the considered levels and cannot be captured by the simple arguments presented here. We will come back 
to this point in section Note however that the distance between £o,cr=o an d £ i,a-=o seem to remain constant as m 



increases which is a manifestation of the presence of the tail correction to the semicircle law (3.2) and is beyond the 
reach of the simplified reasoning we have presented. 

We can next calculate perturbatively the energy of the yrast state in each sector up to the second order in U/A. 
These states can be written as (the singlet and triplet creation operators S\ 1 and T^ 2 (0) have been defined in ( 2.1 2| ) 



and Q2.13Q ) 

l*^ =0) )=4i|0) 

l< =1) > = T+ 2 (0)|0> (3.5) 
Up to the first order their energies are given by 
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£ , a = 2ei + Wlj / 2 
S .i=e 1 +e 2 + V 1 1 ^/2-2XU 

and the second order corrections read (using the constant spacing model for the one-particle levels) 



(3.6) 
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q>/3>2 



(a + /3-2)A 



E 



GO 2 , (^; Q 1 ) 2 (i- ( 5 Q , 2 /2) 



(a-l)A 



(a - 2)A 



^ c o,i 
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a>/3>2 



(a + ,3-3) A 
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a>2 



(^a 2 ) 2 
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1)A (a-2)A 



(3.7) 



Note that any double occupancy in either the initial or the final state, results in a 1/V2 reduction of the transition 
amplitude, hence the factors 1/2 appearing on the right-hand side of the first a nd s econd lines of ([Tt]). These factors 
are however exactly counterbalanced by the IME averages, since one has (see ( 2.14| ) 
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(W^)* = (2U^ + 2U&)* = 8U* 

(O 2 ~~ ~ ~ 



(3.8) 
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rr2,l \2 



AU 2 



The second order contributions for the energies of the lowest levels in each spin sector is therefore given by 
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0,0 
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^ c o,i 
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(3.9) 



The expressions given in equation (|3.9| ) are in very good agre eme nt with numerical data obtained from exact diagonal- 
ization as we show on Figs. || and 0. It is clearly seen from (^^) that the singlet and triplet second order corrections 
differ only by a restriction in the sums which arises in the triplet case because transitions to doubly occupied states 
are not allowed; it is straightforward to show that there are exactly m/2 such transitions. As each contribution in 
second order perturbation theory reduces the energy of the lowest energy state in each sector, these m/2 additional 
transitions will therefore favor a singlet ground-state in the perturbative regime. 

All other transitions give on average the same contribution to £o,o as to £o,i as symmetric and antisymmetric matrix 
elements have the same variance. As the first order corrections do not survive disorder averaging, we can write the 
average energy difference between those two levels in second order perturbation theory as 



T/2 

A 5 w A - 2XU + A— ln(m/2) 



(3.10) 



where A > is a numeric al pr efactor that can be extracted from ( [3.9] ) and the above expression is valid in the large 
m limit. It follows from (|3 .10 ) that in order to align spins, the exchange has to overcome more than just one level 
spacing. Equivalcntly, ( |3.10 ) states th at o ff-diagonal fluctuations increase the energy spacing between the lowest 
energy states of each sector. Equation ( 3.10| ) has been checked numerically and the result is shown in Fig. |[ 

One can also compute perturbatively the splitting induced by the off-diagonal fluctuations between the first a = 
excited state and the a — 1 yrast. As a matter of fact, except for the exchange interaction, all corrections in the 
first two orders in perturbation theory give the same average contributions up to second order contributions which 
exist only for | V I / 1 ^) and correspond to scattering onto a double occupancy. In second order perturbation theory, this 
splitting reads 



£b,i - £ 1,0 = -2AC7 + 



a. a) 



a>3 



2{2a- 3)A 



-2XU 



U 2 

A'— ln(m/2) 



(3.11) 



In particular we see that the splitting induced by the interaction fluctuations favors the spat ially symmetric singlet 
state and opposes the exchange term (A' > 0). Note also that for n = 2, both the splitting (3.11) and the spin gap 







( |3 .10 ) have a similar magnitude . We will see below that this is no longer the case for larger n. Replacing the sum by 



an integral one finds A' w 2 in ( 3.11 ) 



Some remarks are in order here as the case of two particles is somehow special. For n — 2, Uf is a GOE matrix for 
which the number of transitions in each sector is equal to its size. However, as one adds particles, the matrix becomes 
sparser and sparser as the Hilbert space size grows exponentially with the number of particles, whereas the number 
of transitions is a polynomial in n. It is however clear from the perturbative treatment presented above that what 
matters is the number of transitions not the sector size. Generically and for sufficient number of particles, the sector 
with largest number of states has finite (nonzero) magnetization, whereas it is always for a = that one has the 
most transitions and hence the largest probability to find the ground-state. Simultaneously, for increasing number of 
particles, the MBDOS undergoes a crossover to a gaussian shape in the limit n 3> 2LjE£|. It is understood that the 
sparsity of the resulting matrices alone-dpes not invalidate the semicircle law, sparse matrices with uncorrelated matrix 
elements may have a semicircle lawE§c3. However, as noted already, the IME in the TBRIM are highly correlated 
and this apparently drives the MBDOS to the gaussian form. For a very recent and interesting analytical study of 
this crossover, we refer the reader tc£a. Of importance for us is that even for n ^S> 2 one still has a reliable expression 
for the MBDOS in term of n and m that one may use to extract the average energy difference between yrast states 
in the regime of large fluctuations. We will implement this procedure for n > 2 in chapter 

IV. PERTURBATIVE TREATMENT FOR n > 2 

We now discuss the perturbation theory for the yrast states for arbitrary n. These results are of particular interest 
since numerical results for large U are necessarily restricted to small n and one may worry that the large n behavior 
is qualitatively different. In this case, within the perturbative regime, we can show analytically that fluctuations 
reduce the probability of a magnetized ground state for arbitar y n. To estimate the size of U one must consider the 



disorder-ajieraged typical amplitude of fluctuations of the IME (2.4), which has been computed for diffusive metallic 
samplesElreS. In this case the effective static electronic interaction is s tron gly screened and can therefore be well 
approximated by a Hubbard interaction. Then, the variance of the IME (|2.4|) is given by 



« 2 (u2i) -w 2 EE r a (iW^)^)M^)Mj)Mj)^(jW s U) (4.i) 

« i 

In diffusive systems for which l e <C L holds (l e is the elastic mean free path), the wavefunctions can be estimated 
using classical return probabilities as extracted from the diffusion equation and one gets o~(U2'p) = U ~ A/t@. In 
metallic samples the conductance g is very large and even in small quantum dots it is typically of the order of ten. It 
is therefore of interest to start with a perturbative treatment up to second order in the small parameter U/A. Each 
contribution in second order perturbation theory is always negative for each yrast state and we will see, as for the 
case n = 2, that the number of such contributions is larger in the lowest spin sector, thereby favoring the absence of 
magnetization; however additional and more subtle interference effects in the transition matrix elements also appear 
and favor a = 0. Here and if not stated otherwise in the rest of the paper, we will make use of SRS and consider each 
a sector in the a z — block. This means that there are as many particles with up as with down spins, and states with 
different cr's but the same occupancies will differ only in the nature of two-particle bonds between pairs of fcrmions 
on partly occupied orbitals (see Fig. | and the discussion in section ||) . We will also focus most of our discussion on 
the case of an even number of particles, but will eventually generalize our results to an odd number of particles. To 
simplify numerical checks of the perturbation theory we will consider only the case of equidistant one-body orbitals 
e a = (a — 1)A in this section and will discuss generic spectra later on. 

For er 2 = 0, there are an equal number of spin up and spin down fermions and iV(0) = {^n/2^j Slater determinants. 
At U — the ground-state can be written as 

i^>= n ii4j°> ( 42 ) 

a=l,n/2 

Obviously this state has a = 0, as doubly occupied orbitals form a singlet two-particle state. Acting on \F n ) with 
the and T^(0) operators (see ((2.12- 2~13| )) Q and P times respectively 



a, {3 7,5 
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allows to construct a a z — state which is in general a linear combination of Slater determinants of total spin 
a = 1,2, ...Q — 1,Q- One can in pri ncip le represent a complete basis with good quantum numbers a, a z and one- 
particle occupations from the states ( |4.3| ) following the rules : 

• Fermions on the same orbital are singlet paired 

• Fermions on singly occupied orbitals are arbitrarily bonded in pairs, a of the latter being triplet, the rest being 
singlet bonded 

• The triplet bonds combine to maximize the total spin 

While the first rule is imposed by the Pauli principle, the second and third rules-.are a matter of convention. This 
set of rules is similar to the one employed by Kaplan, Papenbrock and Johnsonc3 for the case of n = 4 particles. 
As noted above, the generalization to more particles is not trivial: following the above prescription, one obtains an 
overcomplete basis and one should construct a proper orthogonalization procedure to reduce this basis. In what 
follows however we will compute perturbative corrections up to the second order for only three different states : the 
(7 = and (7 = 1 yrast states (\^q'~ ^} and l^g 17- ^)), and the first a = excited state (|4 , i' T_0 ' > )) For comparison of 
these states the construction of a basis for n = 4 is sufficient. We can write these three states as 

l^ 0) > = \ F n) 



\* i o ) )=Tl /2jn/2+1 (0)\F n _ 2 ) (4.4) 

The difference between the a = 1 yrast state and the first a = excited state lies exclusively in the bond between 
the las t tw o particles : it is a triplet in the first case and a singlet in the second. Up to first order, the energies of the 



states (4.4) are given by 



°0,Q — 


<< 0) 




n ,n 


-1)A 


°i,o — 






~ c 0,0 ' 


-A 


°o,i — 




H\*P) 


— ''O.O ' 


- A - 2XU 



(4.5) 



Without interactions, the latter two levels are degenerate and in first order they are on average splitted only by the 
exchange interaction favoring as usual the spatially antisymmetric triplet state. To calculate the average second order 
corrections, we need to know the number of direct interaction induced transitions which we will call the connectivity 
K and which is calculated in detail in Appendix B. K is a monotonously decreasing function of the total spin and 
in particular the difference between its values at a = and a = 1 is always m/2, independently on the number of 
particles. This decrease of K as a function of a results in a smaller number of second order contributions for states in 
higher a sectors and thus a smaller reduction of the energy of the corresponding yrast state. We will identify below 
the transitions which give the major contributions to the difference in second order shift between the two lowest yrast 
states. The second order correction to the energy of the a = yrast reads 



< } = E E 

a>P>n/2+l 7<<5<n/2 



(VQ 2 (1-<W2)(1- Wg) 
(a + /3 - 7 - 8) A 



+ Y Y 1 a ' fi) rr w A^-n 2 m\n(m) (4.6) 

^ , (a + - 7 - 8 A A 1 ' y ' 

a>0>n/2+l 7<<5<n/2 v ^ ' ' 

Note that the singlet and triplet contributions add incoherently and that trip let t ransition acquire a factor 3 reflecting 



the corresponding number of channels (o~ z = 0, ±1). In order to estimate (^6), the sums can be replaced by four- 
fold integral which gives an homogeneous polynomial of order three in n and m, each term being multiplied by a 
logarithmic correction. In the dilute lim it 1 <C n <C m the m 3 and m 2 n terms drop out exactly and this gives the 
dom inant n 2 m dependence expressed in (E^q). This estimate is also confirmed by numerical evaluation of the sum in 



([16D 



The above formula is found to be in good agreement with numerical data as shown on Fig. |9|. Note that at larger 
number of particles, the dependence of the energies of the yrast states starts to have a linear dependence in U/A much 
earlier, signaling an earlier breakdown of perturbation theory than for small number of particles. We will discuss this 

(2) 

point below. The correction A£q { for the a = 1 yrast can be calculated in the same way and one can show that 
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(2) (2) 

differences between A£p q and A£q { occur first due to denominators differing by ±A as transitions involving the two 
uppermost particles start from the orbitals (n/2, n/2) and (n/2, n/2 + 1) for a = and 1 respectively and second 
due to transitions either increasing or decreasing the number of double occupancies (which only occur for a = 0). 
As noted, the number of such transitions is m/2 and they give a contribution to the spin gap which can be written 
(v = n/ra is the filling factor) 

2A^\n(u) (4.7) 



where A is a numerical factor. This is exactly analogous to the energy difference we found in ( 3.10 ),( [jTl| ) except that 
ln(m/2) has been replaced by ln(^). 

While the term just calculated is easiest to identify, a more important contribution to the spin energy gap comes from 
a more subtle source. There is a certain class of transitions starting from the a = 1 ground-state which have exactly 
the same energy denominator as the corresponding class in a = case (see Fig. [ll]) but the a = 1 transitions have a 
reduced amplitude in comparison to the a = transitions. The corresponding (negative) second order contributions 
will therefore reduce more strongly the energy of the a — ground-state. The a — 1 transitions with this property 
are of the following kind. The a = 1 non-interacting ground state has two partially occupied levels at the top of the 
fermi sea which are triplet-bonded. The relevant transition cause one of these partially occupied states to become 
doubly-occupied while creating a hole in the Fermi sea and a particle above the Fermi sea. For this kind of scattering 
process the number of double occupancies does not change and one can show that the singlet and triplet terms in 
the hamiltonian induce transitions onto the same final state. Correspondingly, the two transition amplitudes must be 
added coherently, and it turns out that this results in a reduced transition probability from 16?7 2 down to 12U 2 (a 
detailed calculation of the amplitude of these transitions is given in Appendix A). The corresponding contribution to 
the spin gap can be estimated as 

A£$ - Afi@ « -4^ / d X dy V — L — « -4^- ln(u) (4.8) 



-0,0 ^ 0il ~ A ,/ - J n/2 "»J-^y + x - s - n/2 

This result is valid in the dilute limit 1 <C n <C m and this contribution dominates the spin gap as soon as the number 
of particles is sufficiently large, i.e. for n > 4. 

As for n — 2 it is straightforward to calculate the splitting induced by off-diagonal fluctuations between the a = 1 
yrast and the first excited a = state : there is no difference in the energy denominators and there is a one-to-one 
correspondence between all second order contributions for these two states, except for the transitions which don't 
exist for a = 1. The latter correspond to scattering onto a double occupancy on the (n/2 + l) th orbital, or from the 
(n/2) th onto a double occupancy on a previously empty orbital. These then are the only contributions to the average 
splitting, which takes the form 



£ °.«=i - £ i.«=° = 2 £ (n+l-2 7 )A + 2 22 (2a-n-l)A "X (ln(m - n)+ln(n)) (49) 

7<n/2 v " a>n/2+l v ' 

and is thus positive. 

We now briefly discuss the case of odd number of particles. The lowest possible magnetization is a = 1/2, and at 
U/A = 0, the yrast corresponds to a singly occupied (n/2 + l) th orbital above a filled Fermi sea (see Fig. [L0| (a)). 
The next magnetization is a — 3/2 and the corresponding U/A — yrast state is represented on Fig. [To] (b). It has 
three single occupancies above the Fermi sea and one of the two bonds between the corresponding particles must be 
a triplet (the choice of the bond is arbitrary). We identified above the dominant second order contributions to the 
spin gap for even n as those which have an amplitude reduction due to partial occupancies in both the initial and 
final states. An example of such a tr ansit ion for a = 3/2 is depicted on Fig. O. From the presented data one sees 



that the expression corresponding to ( |4.8| ) for the case of odd n reads 

tt2 f n/2 ,m/2 ™/2+2 2 

4 X/ dx ^ E -T 2 T«-4^-ln^) (4.10) 

A 7 Jn/2 z> tn / 2 y + Z - X - t 2 A 

and differs from ( |4.8| ) by the boundary values for the sums over z and t. Correspondingly the contribution to the 
spin gap picks up a factor 3/2 and this results in an even odd effect where the gap asymptotically behaves as 
A s « — ABn hi(v)U 2 / A where B = 1 for even and B = 1.5 for odd number of particles. In particular, it is more 
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difficult to magnetize ,a_ system of odd number of ferm ionsE ll wh ich is in agreement with the experimental results 
presented in referencealjS. The above expression ( 4J3 ) and ( [4.10 ) have been checked numerically and the results are 



shown on Fig. [L2j. Both the even-odd dependence and the n-dependence of the gap are confirmed for larger num ber 



of particles n > 3. Note that the processes mentioned above and leading to the scaling expression (4.8) and ( 4.10 ) do 
not exist for n = 2 and 3 in agreement with the data of Fig. [l2[ 

From the second order corrections to the yrast levels in each sector, it is possible to construct an effective Hamil- 
tonian which takes into account the average effect of the off-diagonal fluctuations of interaction. The number and 
strength of second order contributions decreases with increasing magnetization and the relevant contributions are 
those emphasized in this section corresponding to one partially occupied orbital in both initial and final two-particle 
state. For large magnetization a ^> 1 the second order contributions corresponding to the energy difference A^ 
between the lowest spin yrast and the yrast level of spin a can be approximated by the four-dimensional integral 

— j-r — r-r TT 2 [-nil— a i-m/2 r-n/2+a J -J* 

A ( ^4 2 »-4 =4^- dx dy 1 



A Jo Jn/2+a Jn/2-o y + Z — X — t 

^V(? " - ln(i/)) + ^{H~) + \ M2))) (4.11) 
A 2 3 n 3 



We first note that A^ CT ^ is a monotonically increasing function of a. The first term on the right-hand side of (4.11) 
dominates the low-cr behavior. This term is a generalization of the nU 2 / A term giving rise to the spin gap between 
the a = and a = 1 yrasts. Its a 2 parametric dependence results in a delay of the Stoner instability, equivalently in 
a reduction of the strength of the spin-spin exchange coupling 

nil 2 - -. 

— XUS ■ S — > -(XU - A V ——)S ■ S (4.12) 

where A v is a prefactor of order one, weakly depending on the filling factor v. 

For larger magnetization, i.e. w hen t he polarization ratio becomes finite (roughly at a « n/8), starts to be 
dominated by the second term in (4.11) which has a larger a 3 dependence and hence a stronger effect, beyond the 



simple shift of the Stoner instability just mentioned : it results in a saturation of the ground-state magnetization 
for exchange couplings not much stronger than the critical Stoner value. Its er-dependence suggests an higher-order 
effective spin coupling 

Us oc ^-S 3 (4.13) 

which is switched on roughly at a polarization ratio 2a jn > 1/5 (for which \n(2a/n) + 7/31n(2) becomes positive). 
The higher S^-dependence of this effective coupling can also be obtained from a dimensional analysis. The number of 
second order contributions to the ground-state energy in each sector decreases as a 4 for large enough a (see Appendix 
B). When summing over all of these contributions, we must take into account their energy den omin ator, which leads 



to a ~ a rn(<7 ) parametric dependen ce for the second order contributions, in agreement with (4.11). Neglecting the 



logarithmic correction we finally get ( 4.13 ). It is important to note that this latter effective Hamiltonian term is left 
invariant by both SU{2) rotation in spin space and rotation in the one-body Hilbert space. 

The above treatment illustrates the average magnetization decreasing effect of the interaction fluctuations which 
results in a shift of the Stoner threshold to higher exchange strength. Simultaneously, contributions to the fluctuations 
of the ground-state energy around this average in a finite-sized system (quantum dot) can be of the same order of 
magnitude as the average itself, possibly resulting in large fluctuations of the ground-state spin around its average 
value. We therefore close this section with a calculation of the contributions to the variance of the ground-state energy 
arising from the fluctuations of interaction. In first order we get a contribution given by the variance of diagonal 
Hamiltonian matrix elements 

v\(*W\U f \*P)) = 0(n 2 U 2 ) (4.14) 



while the variance of the second order is given by the square of the average contribution (4.6) and is therefore of order 
0(n 4 m 2 £/ 4 /A 2 ). Consequently, these fluctuations are dominated by the second order for U > A/(mn). Of physical 
relevance however are relative fluctuations between ground-states in different spin sectors or at different number of 
particles. Both these quantities ipflpence for instance the distribution of conductance peak spacings for quantum dots 
in the Coulomb blockade regimea'O. In first order, the relative fluctuations between the ground-states with n and 
n + 1 particles are given by 
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(4.15) 

and the relative fluctuations between consecutive (i.e. a and a + 1) yrasts can be written 

- (ZWfc+i 1 - KnZ)^j « VHc/ (4.16) 

where the sums run over occupied orbitals. Both these last two expressions have the same parametric dependence 
on n as they both depend only on the change of one orbital occupancy in the immediate vicinity of the Fermi level. 
In second order, the relative fluctuations between consecutive yrasts can be estimated to have the same order of 



magnitude as the spin gap fl4.8| )-(4.10), and this also gives the contributions to the relative fluctuations between 
ground-states of consecutive number of particles, i.e. 0(nU 2 / A). Consequently, the relative fluctuations will be 
dominated by the first (second) order for U < A/^/n (U > A/y/n). These estimates neglect however the spectral 
fluctuations and are thus valid in the case of a rigid equidistant spectrum only. The variance of the gap distribution 
is however dominated by these spectral fluctuations (which are proportional to the average level spacing A) for both 
Wigner-Dyson and Poisson statistics as long as U < A/y/n. 

V. ASYMPTOTIC REGIME 

In the regime of dominant fluctuations U/A 3> 1, all energy scales (Width of the MBDOS, ground-state energy, 
gaps and splitting between eigenvalues...) become linear in the fluctuation strength U . Most of the properties of the 
Hamiltonian can then be obtained by assuming H ~ Uf, and for random interaction models of this form, the shape 
and width of the MBDOS can be extracted from a computation o£,its variance and higher moments. We begin this 
chapter with a short overview of this method mostly developed ir£3. 

When Uf dominates, Ho and U aV g may introduce a constant shift of the full MBDOS due to the mean field charge- 
charge interaction, a shift of each sector's MBDOS by an amount —XU a(a+ 1) from the mean-field spin-spin exchange 
and a subdominant (0(A/U)) nonhomogeneous modification of the MBDOS due to Ho which is negligible in the limit 
considered here. We thus first consider the MBDOS corresponding to Uf and will introduce later on the only relevant 
mean-field contributions : the cr-dependent shifts due to the exchange interaction. The average shape and width of 
the MBDOS of Uf can be extracted from its moments 

where Uj' J — (I\Uf\J) and |7) refers to a Slater determinant. Taking the average of this expression, we easily see that 
only the even moments of the average MBDOS do not vanish. In the last sum furthermore, only terms with pairs 
of indices occuring twice give a non-zero average contribution, i.e. to compute the moments, one needs to perform 
contractions over the Hamiltonian operators such that 

(I k ,I k+1 ) = (I h I l+1 ) (5.2) 

for a pair of indices (k, I). We can readily calculate the second moment, i.e. the variance of the MBDOS 

M^{a) = ^TrW 2 (a) = (^^-K (a) + (n - l)K x {a) + K 2 {a))W 2 (5.3) 

where we have taken care of the number of matrix elements and different variances of the thre e cl asses of Hamiltonian 
matrix elements mentioned in chapter E3. In the limit of large number of particles, equation (5.3) explicitly expresses 



the dominance of generic two-body IME : since their number is given by K 2(a) ~ n 2 m 2 , their contribution to M^ 2 \a) 
goes parametrically like n 2 m 2 U 2 , whereas the contribution from one-body IME and diagonal matrix elements is 
n 2 mll 2 and n 2 U 2 respectively. This motivates us to neglect the subdominant contributions to MS 2 ' {a) and to use 
the approximation 

M {2) w K(a)U 2 (5.4) 
A calculation of the connectivity K(a) is given in Appendix 
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Higher moments are also easily estimated in the dilute limit 1 <C n <C m. In this case, the contractions (5.2) can be 
performed independently as the probability to create (or destroy) the same fermion on the same orbital is vanishingly 
small. This remains the case as long as the number of creation (destruction) operators in W is smaller than the 
number of particles, i.e. for 2j <C n. A second condition n <C m must also be satisfied, for which creations and 
destructions statistically occur on different orbitals. In this case, higher moments are simply multiple of the second 
moment, with a combinatorial factor reflecting the number of different possible contractions 

M (2j) {a) = {2j - 1)!! [M {2) {a)y (5.5) 

This relation defines a Gaussian MBDOS, and corrections occur only due to higher moments (2j > n), mainly affect 
the tails of the distribution, and vanish in the large n limit. It is remarkable that the order of the moments which fail 
to behave like those of a Gaussian distribution depends almost exclusively on the number of particles, at least as long 
as one restricts oneself to the lowest magnetization blocks away from full polarization. Therefore, corrections affect 
each partial (i.e. er-dependent, or block-) MBDOS in the same way, and we will assume that the relative parametric 
dependence of the bulk of the MBDOS at different a can be extrapolated to the tails. This means that, as in the 
case of a gaussian distribution, the knowledge of its variance full y de termines the MBDOS r-For more details on the 



shape of the MBDOS for random interaction models similar to (2.3) we refer the reader and the more precise, 
very recent treatment given irc3i ■ 

Based on these previous worksliH establishing the quasi-Gaussian shape of the MBDOS, we can now derive a simple 
parametric expression for the average energy difference between yrast states in each magnetization block. Indeed, the 
MBDOS satisfy a scaling law 

E^>E=— (5.6) 

VW) 

which allows to rescale all of them approximately on top of each other. This behavio r is ill ustrated in Fig. [l3] which 
shows both the multiple gaussian structure of the MBDOS and the scaling with y/ K(a) obtained from numerical 
calculations for A = 0. 

The yrast states are distributed in the low energy tail of the partial MBDOS, where the corrections due to higher 
moments are the largest. We have nevertheless seen above that these corrections affect each block's MBDOS in the 
same way (this is true only for not too large magnetizations). Thus the tails undergo the same modification, say for 
a = and a = 1. If we then make the (a priori not justified) assumption that the yrast levels are uncorrelated, i.e. 
that for a given realization of Uf their positions around their respective average value are not correlated, then we can 
conclude that the average distance between two yrast states is parametrically given by the difference of the width 
of the corresponding MBDOS. Assuming, as just discussed, that the tails of the distribution scale with the variance 
with a factor (3 and neglecting contributions arising from Hq, the typical spin gap can be estimated (for A = 0) as 



A? « 0U[y/K(* min ) - y/K(<J min + 1)] (5.7) 

In Fig. Q we show the computed spin gap A^ 7 between the minimally magnetized ground-state and the first spin 
excited level for A = in the limit of dominant interaction, i.e. for Uf. One of the main features emerging from 
the presented numerical data is a strong even-odd effect which is reminiscent of a similar behavior in the limit of 
vanishing interactions. However the origin here is the fluctuating interaction and the energy differences scale as U 
instead of A. As in the perturbative regime discussed in the previous chapter, the occurence of this even-odd effect is 
due to the connectivity K and from Fig. |l4| we see that the probability for a magnetic ground-state is more strongly 
reduced for odd than for even number of particles also in the asymptotic regime. We next note that the gap first 
increases with increasing number of particles before it seems to stabilize above n = 6. We have che cked (dashed and 
dotted-dashed lined in Fig. |lj) that this behavior, which is not captured by the dilute estimate (5/7), is partly due to 



the neglect in (5.3) of nongeneric matrix elements with enhanced variance mentioned above. However, even though 
the exact variance gives a much better estimate, it still underestimates the gap at larger n and we have numerically 
determined that this is due to a strong positive correlation of the ground state energies in adjoining spin blocks which 
is larger at large n. Qualitatively, these correlations are due to the fact that the different block hamiltonians are 
not statistically independent, but are constructed out of the same set of two-body matrix elements. More precisely, 
for a given realization of Uf, all blocks have K{a)N{a) — 0(exp(n), exp(m)) nonzero matrix elements which are 
constructed out of the same set of only 0(m 4 ) different two-body interaction matrix elements. Yrast levels are then 
due to special realizations of the latter inside the blocks. These realizations are presumably not very d iffer ent in blocks 



with consecutive magnetization which results in strong eigenvalues correlations. The above estimate (5.7) which relies 
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only on distribution averages completely neglects these correlations. This is the reason why it underestimates the gap 
at larger n where they are largest. 

The arguments presented in this section are based on estimates for the average yrast energy in each sector extracted 
from the shape and width of the corresponding MBDOS. We have seen in particular that the MBDOS in low spin 
sectors and for a sufficent num ber of particles are almost gaussian with a width given by the square root of the 
corresponding connectivity (B4) \f~K ~ nm. It follows that the ground-state energy in each sector roughly satisfies 
£o.<r ~ nr aU i n the asymptotic regime, whereas in the perturbative regime we found £o,ct ~ n 2 mU 2 ln(m)/A (see 
Equation (4.6)). Neglecting logarithmic corrections (which arise due to the denominators in the second order of per- 
turbation theory) we arrive at the critical border between perturbative and asymptotic regime (radius of convergence 
of the perturbation theory) 



U c ~ A/r 



(5.8) 



Equation (5.8) indicates the breakdown of perturbation theory at a much smaller strength of the fluctuations of 
interaction than previously expected. This is due to the coherent addition of many small second order contributions 
for the perturbation expansion in.the immediate vicinity of the ground-state. A more detailed study of this breakdown 
has been presented in references. 



VI. SPIN POLARIZATION THRESHOLD : DISCREPANCIES FROM STONER'S SCENARIO 

Having established the demagnetizing effect of the off-diagonal fluctuations both in the perturbative and asymptotic 
regimes at A = 0, we now switch on the mean-field spin-spin interaction A > 0. The competition between one- 
body energy, exchange interaction and off-diagonal fluctuations will determine both the average threshold at which 
the ground-state starts to be polarized and the probability of finding a magnetized ground-state at a given set of 
parameters (At// A, U/A). The theory presented in the previous chapters focused essentially on the first aspect and 
we already know that the average threshold for magnetization is increased by non-zero interaction fluctuations. The 
exchange induces energy shifts of —XUa(a + 1) of each sector's MBDOS but has no effect whatsoever on the width 
of the MBDOS. Considering first the asymptotic regime, the average spin gap becomes A s = A^ — XU, where 
A = (5 — (— l)")A/2. In particular, the relative shift betw een the two lowest magnetized blocks is larger for odd 
number of particles, as is the spin gap (See Fig. |l4j). From ( |5.7| ) the average threshold becomes parametrically 

A c ~ V K(a min ) - y/ K(cr min + 1) (6.1) 

From Fig. |A C « 2.5 (3.5) for even (odd) n. Note that as both the spin gap and the exchange are linear in U in 
the asymptotic regi me, t his average threshold is [/-independent. This is no longer the case in the perturbative regime. 



As shown in section IV, the perturbative spin gap can be approximated by A S (U) — A ~ BnU /A, where we recall 



that B = 1 (1.5) for even (odd) n. We then get 

\ nit 



Ac((7) - (Oo ~ Bn U/A (6.2) 



where we used the critical (Stoner) exchange strength {U^'p)o = A/2. Once this threshold is reached, the spins start 
to align, but in contrast to the Stoner scenario, full polarization is not achieved at once, because of a parametric 
decrease of the second-order contributions from off-diagonal fluctuations as a is increased. From the perturbative 
treatment presented in section IV a S 3 term takes over at large spin which induces saturation of the ground-state 



spin. The mechanism for the appearance of that term is a reduction of the probability for transitions from or onto 
partially occupied orbitals with respect to transitions from doubly occupied orbitals onto empty orbitals. Off-diagonal 
fluctuations result in two effective Hamiltonian terms ~ S ■ S and ~ S 3 and the second term influences the system's 
magnetization properties at large spin, but before full polarization. Neglecting logarithmic corrections in n, m and a 
and for a given A = X C (U) + SX (i.e. A measures the distance to the Stoner threshold), A and U, the magnetization 
will saturate at a value 

CTmax ~ SX— (6.3) 

This is a major modification of the Stoner scenario for which once the magnetization threshold is reached, full 
polarization of the electrons is achieved at once. The presence of off-diagonal fluctuations, no matter how weak, induces 
this saturation, as their relative weakness will eventually be counterbalanced by the larger parametric dependence in 
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a of the number of second order contributions at large a. We stress that this saturation is entirely induced by the 
off-diagonal fluctuations and does not depend on any modification of the one-body density of states at larger spin. 

We next show on Fig. |l5] the behavior of the spin gap between the two lowest yrasts as a function of U/A and 
for different values of A. The variance-of the gap distribution is of course unaffected by the exchange and we already 
know that the probability P(a > 0]B of finding a magnetized ground-state is reduced by the off-diagonal matrix 
elements. This probability will eventually saturate above a finite value of U/A, since the width of the gap distribution 
is proportional to its average ~ U / AE2l. This is shown on Fig. |l5| where the error bars reflect the width of the gap 
distribution. Their linear increase with U means that the fraction of negative "gaps" (contributing to the probability 
of being magnetized) is constant with U. The same behavior is characteristic of gaps between higher consecutive 
yrast, which results in a [/-independent behavior of P(a > 0) at large U. 

Finally, P(a > 0) is shown on Fig. [iTJ as a function of the exchange strength XU/A for different values of U/A and 
different distributions of one-particle orbitals. This figure shows a clear demagnetizing effect of the fluctuations of 
interaction except below the Stoner threshold in the case of an equidistant spectrum. We recall that the demagnetizing 
effect is in fact only an average effect, and that for an equidistant spectrum, U may for particular realizations reduce 
the level density at the Fermi level, thereby favoring the appearance of a higher spin ground-state as can be seen on 
Fig. [if] for small interaction fluctuations U/A — 0.1 and small exchange strength XU /A < 0.5. In the two other cases 
of a randomly distributed and Wigner-Dyson one-body spectrum, fluctuations of interaction always reduce P(a > 0). 
At larger U, the dependence on the orbital distribution is rather weak, as shown in Fig. [l7| Note in Fig. [TJ the 
bending of P(a > 0) above the onset of magnetization which is a clear difference from the Stoner behavior : even 
at quite large exchange, P(<r > 0) remains smaller than one. From these data, we define an average magnetization 
threshold X C (U) for which P(a > 0) = 0.01 and extract from Fig. [TJ| the additional exchange strength <5A necessary to 
achieve P(a > 0) = 0.5. The results are shown on Fig. [l^ and indicate a linear increase of SX with U which illustrate 
the demagnetizing effect of the off-diagonal fluctuations : a stronger exchange than predicted by a simple Stoner 
picture is necessary to have even a weak nonzero ground-state magnetization probability (see Fig. |l6| ), moreover an 
even stronger one is necessary to achieve a significant probability. All this is in qualitative agreement with equation 
(5.3). A direct numerical check of this equation would however require a much larger number of particles, beyond 
today's numerical capabilities. 



VII. REAL SPACE MODELS 



It is now evident from the results presented above that fluctuations of IME introduce a new energy scale. In 
addition to the Stoner parameter XU/A, the ratio A between the exchange and the interaction fluctuations gives 
a second relevant parameter for the emergence of a ferromagnetic phase. We therefore turn our attention to the 
microscopic computation of the magnetization parameter A for standard solid-state models. This will allow us to 
estimate the strength of the demagnetizing effect of off-diagonal fluctuations in more realistic situations. We consider 
Anderson lattices whose one-body Hamiltonian is given by 

Here restricts the sum to nearest neighbors, and Wi € [-W/2; W/2] where W is the disorder strength. We study 
interaction potentials of the form 

U{% - j) = U 5(i - j)+Ux/ \n - fj\ (7.2) 

i.e. for IA\ — we have a pure Hubbard interaction whereas IA\ ^ gives a long-range interaction. Microscopically, A 
is given by the ratio of the average exchange term 

(O = E U ( l - 3)r a {i)%{j)M)Mi) (7-3) 



and the r.m.s. of the distribution of IME (|2.4|) . By definition the average in equation (7.3) is performed over 
wavefunctions close to the Fermi level. Fig. |19| and |2(] show the disorder dependence of A, for a pure Hubbard 
interaction U\ = on two- and three-dimensional lattices respectively and for different linear system sizes. The data 
have been obtained from averages over 30 wavefunctions in the middle of the Anderson band E = and for 10 (L = 80 
in 2D and L = 15 in 3D) to 200 (L = 10 in 2D and L — 6 in 3D) disorder realizations. In both dimensionalities we 
can distinguish three regimes : (I) At low disorder, the one-electron dynamics undergoes a crossover from ballistic 
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to diffusive regime as the linear system size is increased beyond the elastic mean-free path l e ~ (V/W) 2 . In the 
ballistic regime l e 3> L, wavefunctions are plane- waves. In this case, a Hubbard interaction gives A ~ L 2 , since the 



RMS 



TP' 5 



> and (DJJ) ~ L-i, whereas once the diffusive regime is reached, one expects A ~ A/(A/g) ~ 30. 



In the crossover between these two regimes, contributions from gaussian modes (those corresponding to \i — j\ < l e 
in (4.1)) may dominate the fluctuations of the IME but eventually vanish as one increases L as they are weighted 
by a factor {l e /L) c \^. Presumably these contributions still affect our data in region (I). (//) In the regime of 
intermediate disorder, both off-diagonal fluctuations and exchange are increased by disorder, and apparently they 
compensate each other, resulting in a L-independent A w 4, in 2D. We expect that this behavior will hold as one 
further increase the system size. We indeed numerically estimated the elastic mean free path at W/V = 5 from the 
distribution of inverse participation raticH and found a value l e 4. The gaussian modes are thus weighted by a 
prefactor 1/400 for L = 80 and have therefore only a marginal influence on the fluctuations of the IME, so that one 
may reasonably assume that finite-size effects have only a marginal influence on the data presented in Fig. [1^ in region 
(II). In the three-dimensional case, it even seems that A decreases as the system size increases in the intermediate 
regime W/V <G [8, 17], however this is due to the quite small linear system sizes considered here, and once one reaches 
L 3> Z e , A should saturate at a finite, but quite small value. It is interesting to note that the upper border of this 
intermediate regime is quite close to the critical disorder value for the Anderson localization transition. (Ill) In the 
regime of strong disorder, one-particle wavefunctions are strongly localized on fewer and fewer sites, the off-diagonal 
fluctuations are sharply reduced (due to quasi selection rules discussed in chapter and again exchange dominates. 
Note that eventually, the latter disappears also, but at a lower rate than the fluctuations. These results indicate that 
at an intermediate disorder strength, off-diagonal fluctuations may be strong enough to play an important role for 
the magnetization properties of the ground-state. 

We next evaluate the influence of the long-range part of the interaction. The average exchange interaction (7.3) 
term is given by an average taken over one-particle wavefunctions close to the Fermi level. Due to their orthogonality, 
taking this average over the full set of wavefunctions gives a 5- function and only on-site contributions. This averaging 
procedure is however only justified if the one-body dynamics is described by Random Matrix Theory (RMT) for 
which the structure of the eigenstates is homogeneous all through the spectrum. RMT however describes systems 
which are of interest here only inside an energy window given by the Thouless energy around the Fermi leveled so 
that the average over wavefunctions close to the Fermi level leads only to a more or less sharply peaked function of 
(fi — fj). There are also contributions to the exchange from the long-range terms, but still we expect that the average 
damps them with respect to their contribution to off-diagonal fluctuations (This damping of course depends on the 
disorder strength.) which are of the same order of magnitude as the short-range contribution up to distances of the 
order of l e c3. This means that we expect a decrease of A upon increase of the interaction range. The validity of this 
reasoning is illustrated for the two-dimensional case on Fig. ^l] where we plot the evolution of A for different disorders 
as the long-range part of the interaction becomes more and more important. Clearly, A decreases as the range of the 
electron-electron interaction increases, and therefore the Hubbard results presented on Fig. |l^ and ^ give an upper 
bound for A. One thus expects the demagnetizing effect described in this paper to be more efficient at low filling when 
the screening length exceeds the elastic mean free path. 

In finite-sized systems like quantum dots where impurity scattering is weak but wavefunction fluctuations are 
induced by chaotic scattering at an irregular confining potential, standard estimates give A rs g for a short-range 
interaction, whereas in the (unphysical) limit of an infinite range interaction U(r — f 1 ) = U one gets A = lo. 
Therefore, and as g is not too large in such systems, it is a priori not justified to neglect the effect of off-diagonal 
fluctuations, as they should at least strongly suppress the probability of finding ground states of larger spin beyond 
few (ss 2, 3) polarized electrons. It has even been proposed by Blanter, Mirlin and Muzykantskii that in confined 
systems the accumulation of charge at the surface of confinement leads to stronger fluctuations of screened Coulomb 
interaction matrix elements ~ A/y^ which would give A ~ y/g. As in quantum dots g is of the order of up to few tens, 
this would bring A down to values where the demagnetizing effect of off-diagonal fluctuations plays an important role. 
All this illustrates the relevance of off-diagonal fluctuations for the magnetization properties of the ground-state in 
regimes of intermediate disorder and for poor screening of the electronic interactions - presumably, for low electronic 
densities p for which the distance between electron is smaller than the elastic mean free path p 1 / < l e . 

Assuming still U ~ A/g, the shift of the Stoner threshold is quite small, of the order 0(A/g). This is so, as the 
model we consider is valid only in an energy window of the order of the Thouless energy E c = gA aroun d the Fermi 
level, so that it is quite natural to set n,m ~ g. At larger magnetization however, the second term in (4.11) takes 
over and induces a significant reduction of the ground-state spin when the latter becomes comparable to g with a 
prefactor depending on the strength of the average exchange. This term strongly modifies the Stoner scenario as it 
induces magnetization saturation above the magnetization threshold and full polarization can be achieved only once 
a second, significantly larger, threshold is reached. 

We finally describe an experimental setup that allows to gain important information on the ground-state spin 
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of two-dimensional lateral quantum dots in the Coulomb blockade regime. The experiment was proposed irfi and 
consists in applying an external magnetic field in the plane of a lateral, two-dimensional quantum dot. Because of the 
two-dimensional nature of-the dynamics, we assume that an in-plane field has no orbital effect so that it introduces 
only a Zeemann couplings. The difference in ground-state spins for consecutive number of electrons can then be 
determined experimentally by studying the motion of Coulomb blockade conductance peaks at very low temperature 
T < A as the in-plane magnetic field is increased. The resonant gate voltage is given by a difference of two many- 
body ground-state energies eV™ — £^+1 ~ an d it is always the difference of an even-odd pair (i.e. of ground state 
energies corresponding to one even and one odd number of electrons on the dot). Upon application of an in-plane 
field, the peak position behaves like 

eV g n (B) = -£°+ gl x B Bba z {n) (7.4) 

Sa z (n) is the magnetization difference between the two consecutive ground-states which can therefore be extracted 
experimentally from the motion of conductance peaks in an in-plane field. At minimal magnetization one has a 
sequence of ground state spins a z — 0, 1/2, 0, 1/2, 0, 1/2... (For odd n, and due to SRS, the a z = 1/2 and a z = —1/2 
ground-states are degenerate so that an arbitrarily weak in-plane field aligns the spins and one always has a z = 1/2), 
therefore 5a z (n) = (— 1)"/2 and one has \dV g /dB\ — g^is/^. As B is increased the ground-state with even number 



of electrons is most likely to magnetize first (this is because of both the even-odd effect mentioned in section IV 
and the larger kinetic energy cost to flip one spin for odd number of electrons), exactly reversing the slope of two 
consecutive peaks; then as the field increases further the odd state will likely flip, restoring the original slope. As long 
as consecutive ground states never differ by more than one unit of spin the absolute value of the slope will remain 
constant as the system polarizes. Correspondingly, if all slopes are constant \8V g /dB\ — it is very likely that 

no ground-state is magnetized (this would for instance require a sequence <r z = 1, 1/2, 1, 1/2, 1, 1/2... We do not 
see any obvious reason why all even n ground-states should have a z — 1 while at the same time none of the odd n 
ground-states are magnetized). 

However if there exist many magnetized ground states, then the probability to find pairs of consecutive ground- 
states with a larger difference in magnetization |&r z (n)| > 1 increases and one expects a range-of slopes to occur. 
Then the corresponding peak heights may be strongly reduced by the spin blockade mechanisnH, which should be 
easily visible experimentally. This argument ne glec ts changes in the g-factor of the electron with changing n, which 
presumably are slow. This is illustrated on Fig. £2] where the peak positions are drawn as a function of the Zeemann 
coupling for A = 1 and 5. It is clearly seen that at weak A, \dV g /dB\ is constant and corresponds to a minimal Sa z , 
while a larger A gives different slopes in agreement with the above reasoning. Note also in this latter case, that peaks 
evolve in parallel at weak magnetic field, indicating the successive addition of two spins oriented in the same direction. 
This feature is absent of the weak-exchange (right-hand) side of the graph for which the ground-states arc obtained 
by piling up electrons on the orbitals according to the Pauli prescription. This results in a minimal ground-state 
spin, a sequence Sa z (n) = (— 1)"/2 of magnetization differences between consecutive ground-states and a motion of 
neighboring peaks in opposite direction at low field. 

Recently, experiments in this direction have indeed been performed, which have given-setious evidences for the 
occurence of partially (but weakly) magnetized ground-states with few polarized electronsEj'E3. This means that the 
off-diagonal fluctuations are not dominant, in agreement with the above considerations giving a large ratio A ~ g 
between the strength of the average exchange and the off-diagonal fluctuations. Therefore the perturbative treatment 
presented in section 



(large dots) or in the 



IV is expected to correctly describe semiconductor quantum dots cither in the diffusive regime 



ballistic regime with chaotic boundary scattering (smaller dots). We note in this respect that from 



recent experiments on the distribution of conductance peak spacings, the parameters A and U have been tentatively 
extracted and, seem to indicate a conductance g ~ 6 for which off-diagonal fluctuations should give a non-negligible 
contribution^. Finally, we note another interesting experimental result which is the apparent absence of suppression 
of the conductance peak found in some cases for larger spin difference between consecutive ground-statespl. This is 
in major disagreement with the spin-blockade mechanism proposed by Weinmann, Hausler and KrameiEJ and is yet 
to be understood. 



VIII. CONCLUSIONS 



In this article we have illustrated how fluctuations of the interaction matrix elements tend to reduce the ground- 
state magnetization, both when they can be treated perturbatively (the regime which is relevant for condensed matter 
physics) and in the asymptotic regime where they give the dominant terms in the Hamiltonian (which is relevant 
for nuclear physics). The mechanism behind this effect is in a way similar to the Stoner picture where itinerant 
ferromagnetism occurs due to a larger number of diagonal interactions at low magnetization. As in a mean-field or 
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self-consistent approach, each of these terms gives a positive contribution (for a repulsive interaction), this directly 
favors spin polarization. Similarly, we have shown that interactions induce more transitions in the low spin sectors. 
Each of these transitions gives one contribution in second order perturbation theory which this time is however negative 
(both for attractive and repulsive interactions) if one considers the lowest level in each sector, and this therefore favors 
a low spin ground-state. In the perturbative regime, we have seen that these fluctuations induce two terms in an 
effective Hamiltonian formalism : a S ■ S term which simply induces a small shift of the Stoner threshold and a 
second l^^-term which is switched on at larger magnetization where it results in a saturation of a. This is a major 
qualitative modification of the Stoner scenario : even neglecting discrepancies in the one-body density of states, full 
polarization is not achieved once ground-state magnetization has been triggered by the exchange interaction. The 
latter must indeed also overcome the ~ IS 1 ! 3 term, which requires an even larger exchange. The strongest effect of 
off-diagonal fluctuations occurs in the large spin regime a = 0(g) where a mean-field picture overestimates the value 
of the ground-state spin. 

From the point of view of nuclear physics, our analysis of the regime of large fluctuations, based on a study of the 
many-body density of states, clearly indicated a strong bias toward low angular momentum ground-state. We have 
not explained however why numericaLresults indicate an almost 100% predominance of o~ = ground-state for models 
similar to the one we have studiedBEJ, and this question is still open. 

Our findings should finally be put in perspective with the Renormalization Group (RG) treatment for disordered 
interacting electronic systems of FinkelsteintEI. In his approach, one indeed finds that the RG flow renormalizes the 
ferromagnetic spin-spin coupling to larger and larger values, possibly indicating the occurence of a ferromagnetic 
phase due to the combined effect of disorder and interaction. The perturbative treatment we presented in section 
IV did not allow us to find any contribution favoring a higher spin and this apparent disagreement between the RG 
approaches and ours is at present not understood. We note however that it has been suggested that the divergence 
of the exchange coupling induced by the RG flow could indicate a crossover to the singlet-only universality classlElEj. 
In this respect it is worth noticing that the scattering processes in the singlet and triplet channels as defined in the 
present work have coupling constants corresponding to the sum and the difference of the couplings T and T 2 as defined 
ir£2l respectively. It can be checked that the ratio a(W)/a(V) = (T + T 2 )/(T — T 2 ) satisfies the same RG equation as 
the exchange coupling (72 = T 2 /z ino) so that the triplet channel vanishes at the same rate as the (ferromagnetic) 
exchange flows to strong coupling, which may indicate a cancellation of the ferromagnetic instability by the effect 
studied in the present paper. 
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APPENDIX A: 



Under a rotatio n in spin space, the triplet operators in equation ( 2.15 ) are brought into one another, whereas the 
singlet operators (2.12) are left invariant. SRS on the other hand implies a number of interaction-induced two-body 
transitions which is invariant under such a rotation. SRS can be easily check ed fo r initial states without double 
occupancy, and it is equally easy to convince oneself that the singlet operators (2.12) is spin conservi ng. F or initial 
states with double occupancies however, the triplet operators (first three terms between brackets in ( 2.10] )) are not 
individually SRS but must be considered as one single spin conserving triplet operator. To check this one ac ts on a 
four-particle state with two double occupancies (which has thus a = a z = 0) with the full triplet operator of ( [2.10 ) 



1 - .1' -sT c <L c <5,tI°) 
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(Al) 



This is the sum of two products of two singlets, and it is obviously spin conserving. Moreover, acting on the same 
initial state with a singlet interaction operator gives 



The two final states in ( |Al| ) 



2~( c i,t c ki ~ c l,i c /3,T)( c 7,i c <5,T ~ c 7,T c <u) C 7,i C 7,t c il c !Tl°) 
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t J \(J J 



t j 

i, 



(A3) 



are not orthogonal to each other and it can easily be checked that the sum 2|#) + |*'> = V3|$> is a singlet. |<i>) is 
normalized and the factor \/3 appears in the case of an initial state with double occupancies for which the n umb er of 
triplet transitions is reduced by a factor three (without double occupancies, the three triplet operators in ( 2.10 ) are 
individually SRS). Thus the total transition probability is kept constant. Finally, |<I>) is orthogonal to |\&'), i.e. singlet 
and triplet channels give transitions onto orthogonal states. Their contribution to second order perturbation theory 
will therefore add incoherently and give a transition probability 



z(v^r+(wz : ;r = i6u 



(A4) 



A calculation going along similar lines shows that if one of the final orbitals (e.g. a or f3) is partially occupied, the 
transition probability is reduced by a factor 1/2. 

We also calculate the transition probability for one singly, one doubly occupied initial orbitals. The initial state to 
consider is 



i*±) = ^(<T c L ±c L c kT) c L c i T io> 



(A5) 



The label ± refers to either a triplet or a singlet a z — two-particle state on the orbitals a and f3. The transition 
amplitudes can be calculated in the same way for both cases, and we restrict ourselves below to the triplet case with 
Note that this latter state can be brought via a rotation in spin space onto the following <r z = 1 state 



I*, 



C a,T C /3,T C <5a 



t J 



and that the calculations to be presented below give the same transition amplitude for both 
thus fully SRS. Acting on |^+) with a singlet interaction operator gives 



(A6) 



and l^,) and are 



c /3,T c 5,i)l*- 



(A7) 



which once again is SRS. Two remarks are in order here. First, the above transition amplitude has picked up a factor 
1/V2 with respect to the case where the initial state has two double occupancies. This is due to the vanishing of 
one transition "channel", as the (3 orbital is only singly occupied, and will result in a factor 1/2 for the transition 
amplitude. Note that this factor is counterbalanced by a twice larger number of transitions for the case considered 
here, as one has the freedom to destroy (or create) a particle on the a th or the /3 th orbital. Secondly, doing the same 
calculation with a triplet operator acting on the singlet initial state 1^+) is not SRS per se, but once again requires 
to consider the a z = ±1 triplet operators, as we did above for the case of two double occupancies. 

Only in the situation where both initial and final states correspond to partially occupied orbitals does one get an 
uncompensated reduction of the transition amplitude with respect to the above case of doubly occupied initial and 
empty final orbitals. As we are now going to show, this results from the coherent addition of the triplet and singlet 
transitions which lead to the same final state. The initial state is e.g. 
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V2 

and one acts on it with the operator 

1 



0± = ^(4,T c m4 ± C 7,1 c Lt)( c /3,1 c 5,T ± c /3,T c 5,l) (A9) 
A straightforward calculation gives the same result for both singlet and triplet operator 

o±\* in ) = ^(<t c L + c \,i c hHA,i\°) ( A1 °) 

The result is SRS, i.e. gives a four-particle final state with a = 1, cr z = 0. The key point here is that both singlet 
and triplet channels go to the same final state. Thus one gets the corresponding second order transition probability 
by adding their amplitude coherently. The average transition probability reads then 

\(v™ + w%*y = (ag + u s ^y = 2u 2 (ah) 



The same contribution arises from the interchange [3 <-> /i in the operator (A9). The above transition probability 
comes therefore with a factor two. This is because once the particles are triplet (or singlet) paired on different orbitals, 
transitions become distinguishable. In addition, one has to consider the four a = 1 transitions induced by (s =f , J.) 

O s = 4,s C Ls C /3^C5, s Or C^C^C^sC^s (A12) 

which together give a transition probability 2V 2 = 8U 2 . The total transition probability is then 12U 2 instead of 16U 2 



in the case of doubly occupied initial and empty final orbitals (A4). The same result is obtained in the case of a 
singlet initial state 

1**0 = -j^IAa - iA,MAi\°) (Ai3) 



for which singlet and tripl et transitions (A9) and (A10) also add coherently, resulting in the same reduction of the 
transition amplitude (All). Consequently, and for n > 3, the splitting between the a = 1 yrast and the first a = 
excited state is much smaller than the spin gap between the two lowest yrasts. 

Note that these calculations must be modified in realistic systems for which V 2 — W 2 does not necessary hold. For a 
purely local (Hubbard) interaction with time-reversal symetry for instance, one has V 2 — ^ W 2 as the antisymetrized 
matrix elements vanish exactly. Then, the reduction in transition probability occurs due to the vanishing of singlet 
transitions as one goes to larger magnetizations. 

APPENDIX B: 

We first calculate each sector's connectivity K which is the number of basis states directly connected to an arbitrary 
initial many-body state of a given sector, alternatively the number of nonzero matrix elements per row (or column) 
of the Hamiltonian matrix. We saw in Appendix A that some transitions have increased weights, in particular triplet 
transitions involving two doubly occupied orbitals pick up a factor y3 that is absent of all other transitions. In 
absence of double occupancies however, these transitions are replaced by three times as many triplet transitions so 
that the total transition probability is conserved. The latter quantity is in fact the physically relevant one as it 
appears in second order perturbation theory and determines the scaling of the MBDOS in the regime of dominant 
fluctuations. We therefore calculate the weighted connectivity, where the number of transitions are multiplied by the 
square of their relative amplitude. With this defi nition and for the case we are considering of spin 1/2 particles, the 
connectivity is constant within one sector. From ( pTOD , K is the sum of a singlet and a triplet channel contribution 



which differ only in that the former allows transition from and to double occupancies. For a — we may consider 
the U = ground-state as our initial state. It is easily seen then that the number of directly connected states can be 
expressed as a sum over four contributions K = K + K\ + K s + K t which for a — are given by 
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The first term corresponds to trivial diagonal transitions, the second one to partially off-diagonal transitions changing 
a single one-body occupation, while the third and fourth one s cor respond to generic two-body transitions induced 



respectively by the singlet and triplet interaction operators in (2.1C). Note the discrepancy in the prefactor 3 between 
K s and K t due to the enhancement of triplet transition amplitude discussed in Appendix A. As a increases, some 
singlet transitions are replaced by additional triplet transitions but some other disappear which we are going to 
identify. The connectivity at full polarization is also easily calculated as there are no long singlet transitions and 
particles may be considered spinless. One has K = Kq + K\ + Kt 

2 - n) (B2) 
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= n/2) 
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= n/2) 


= n 
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= n 



The connectivity difference between minimal and maximal polarization is therefore given by 

3 o 4 

77 777 l77 

K(a = 0) - K(a = n/2) = — + 0(n 3 , m 3 , n 2 m, nm 2 ) (B3) 

8 16 

For a finite magnetization, we may represent the lowest level in the sector as a sea of n/2 — a double occupancies 
separated from m/ 2 — n/2 — a empty levels by a layer of 2a singly occupied levels as depicted on Fig.E3| The transitions 



included in (Bl) and that are now forbidden correspond to singlet transitions involving either initial or final states 
with at least one scattering particle in the 2a layer. These transitions can be classified as : 

1) transitions from a double occupancy in the Fermi sea onto the 2cr-layer (Fig.|24| (a)) 

2) transitions from the 2<r-layer onto a double occupancy in one of the m/2 — n/2 — a empty levels (Fig.|24| (b)) 

3) One- and two-body transitions within the 2c-layer (Fig.^J (c)) and 

4) Two-body transitions from the 2cr-layer, one of the particles being transferred to a new orbital in the 2cr-layer, the 
other one onto one of the m/2 — n/2 — a empty levels (Fig.p4| (d)) 

A simple counting of the number of these transitions finally gives 

K{a) = K(0) - [a(2a - l)(n/2 - a) + a(2a - l)(m/2 - n/2 - a) 

+ (2(7(2(7 - 1) + a(a ~ l) 2 {2a - 3)/2) + a(2a - l)(a - l)(m/2 - n/2 - a)} (B4) 



which in particular correctly reproduces the difference (B3). It is easily checked (e.g. from (B3) that the ratio 



K(a — n/2)/K(0) = 1 — Av + Bv 2 is a function of the filling factor v — n/m only. Also it is remarkable that the 
connectivity difference between a = and a = 1 is m/2 for any number of particles. 
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FIG. 1. Distribution of off-diagonal interaction matrix elements (2.4) for a one-dimensional lattice model with nearest and 
next-nearest neighbor hopping and a Hubbard interaction (Details of the model can be found in referenceEj) and for different 
strength W/V of the disordered potential W/V = (thin solid line), 1 (thick solid line), 2 (dashed line) and 3 (dotted-dashed 
line) . 
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FIG. 2. Representation of <r z = many-body states with a = 2 (a) and (b). These two states differ only by the nature of 
the two-particles bonds connecting pairs of fermions on partly filled orbitals which are either singlets or triplets (the nature of 
the bonds is indicated by the letter s or t). As fermions on doubly occupied orbitals can only be singlet paired, they cannot 
provide for a nonzero spin. Together with SRS, this forbids the scattering from a triplet bond configuration onto a double 
occupancy, so that the rightmost state (c) can only be coupled to the a = state (b). 
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FIG. 3. Density of states for n — 2, m — 12 a — (empt y cir cles) and a = 1 (full circles) computed from 5000 realizations 
of 14 f. The full lines give the corresponding semicircle law (3.2). Tails develop due to the finiteness of the Hilbert space size 
(N(Q) = 78 and N(l) = 66). 
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FIG. 4. Average number p of levels in the spin gap between the two yrasts for n — 2, divided by the total number of levels 
N(0) = m(m + l)/2 in the a = sector as a function of the number m of one-particle orbitals. The dashed line shows the 
dependence p/N ~ 1/N in agreement with a m-independent number of levels in the gap. 
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FIG. 5. Average spin gap between the two lowest yrast levels (circles) and average splitting between the first excited a — 
level and the a = 1 yrast level (triangles) of W/ for n = 2 as a function of the number of one-particle orbitals m. The data 
show almost no m-dependence in agreement with (3.3). 



O 

o 
o 
CO 
I 

g 

o 
d 




-2.0 
0.00 



0.05 



0.10 0.15 

U/A 



0.20 



0.25 



FIG. 6. Ground-state energy for the Hamiltonian ( pi.S| ) with A = at n = 2, a — 0, m = 12 (circles) and m = 16 
(squares) as a function of the strength of off-diagonal fluctuations U/A. The solid lines indicate the perturbative result 
{eo,o(U) — eo,o(0))/A = A(U/A) 2 with a numerical coefficient determined by (BJJ) A = —21.12 and —27.56 respectively. 
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FIG. 7. Ground-state energy for the Hamiltonian (2.9) with n = 2, a = 1, m = 12 (circles) and m = 16 (squares) as 
a function of the strength of off-diagonal fluctuations U/A in absence of exchange interaction. The solid lines indicate the 
perturbative result (eo,i(U) — eo,i(0))/A = A(U/A) 2 with a numerical coefficient determined by (3.9)A = —17.5 and —23.91 
respectively. 
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FIG. 8. Spin gap between the two lowest yrast states for n = 2 and m — 16 as a function of the strength of off-diagonal 
fluctuations U/A. The solid line give s the perturbative result from (3.1C), giving (eo,i{U) - e ,o(t/))/A = 1 + A(U/A) 2 with a 
numerical coefficient determined by (3.10), A = 3.66. 
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FIG. 9. Energy of the a — yrast state for m = 16, 
function of the strength of off-diagonal fluctuations U/A. 

giving (eo,a(U) - eo,o(0))/A = A(U / A) 2 with a numerical coefficient determined by (|4.6|)A = -27.56 (n = 2), -144.75 (n = 4), 



n — 2 (circles), 4 (squares), 6 (diamonds) and 8 (triangles), as a 
The solid lines give the pert urbative results extracted from (4.6), 



-281.09 (n = 6) and -373.57 (n : 
the large U/A linear regime. 



Note that the breakdown of the perturbative expression coincides with the emergence of 
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FIG. 10. The two lowest yrast states a = 1/2 (a) and 3/2 (b) for odd number of fermions. 




FIG. 11. Left: Transitions that have a different transition amplitude for a — and a = 1 and thereby give the dominant 
contribution to the spin gap between the two lowest yrast levels for even number of fermions. Right: Corresponding transitions 
for odd number of fermions giving different transition amplitudes for a = 1/2 and a = 3/2. 
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FIG. 12. Rescaled spin gap between the two lowest yrast states for m = 16 and n — 2 (circles), 3 (squares), 4 (diamonds), 
5 (triangles up), 6 (triangles left), 7 (triangles down) and 8 (triangles right). Symbols corresponding to odd (even) n are filled 
(empty). The scaling parameter satisfies a — 1 (1.5) for even (odd) number of particles (see text). The scaling holds quite well 
already for a small number of particles n > 4 and the even-odd dependence of the gap confirms the theory presented in the 
text. Inset : spin gap before rescaling for the same cases as above. Lines corresponding to odd (even) n are dashed (full) to 
stress the even-odd dependence of the gap. Note that for larger n, the gap starts to have a linear dependence above U/A £s 0.1 
indicating the border between perturbative and asymptotic regimes. 
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FIG. 13. Density of states for the Hamiltonian Uf with n = 6 particles and m = 16 orbitals, corresponding to the magneti- 
zation blocks a = — 3 (solid line), -2 (dotted line), -1 (dashed line) and (dotted-dashed line). Inset: rescaled density of states 
showing the approximate scaling in E/K l ^ 2 U . 
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FIG. 14. Dependence of the finite-size spin gap in the number n of particles in the regime of dominating fluctuations l//A> 1 
and A = 0. Points correspond to numerical results for m = 10 (full circles), 12 (empty squares), 14 (full diamonds) and 16 
(empty triangles). For the case m = 16 and 1000 Hamiltonian realizations, the error bars indicate the r.m.s. of the gap 
distribution while the dashed and dotted-dashed lines show the numerically computed variances (Left-hand side of eq. (5.4)) 
for the full Hamiltonian and after setting to zero non-generic interaction matrix elements respectively. 
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FIG. 15. Evolution of the energetic distance between the two lowest yrast for n — 5 and m — 12 as a function of the strength 
of interaction fluctuations U/A. Shown are curves corresponding to an exchange XU/A = 0. (circles), 2 (squares), 4 (diamonds) 
and 6 (triangles). The error bar indicate the r.m.s of the distribution of A s /A. 




FIG. 16. Probability for a magnetized ground-state as a function of the exchange XU/A for 5000 realizations of Hamiltonian 
( ^.9|) with n — 6, m — 12. Three cases with equidistant (a), GOE (b) and random (c) one-body spectra are shown. Different 
curves correspond to different fluctuations of IME : U/A = (solid line), 0.1, 0.2, 0.3. ..0.8, (symbols, from left to right). 
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FIG. 17. Comparison of the ground-state magnetization probability for n = 6 and m = 12 at (7/ A = (lines) and 0.5 
(symbols) for equidistant (solid line and circles), Poissonian (dashed line and squares) and Wigner-Dyson (dotted-dashed line 
and diamonds) orbitals distribution. 
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FIG. 18. Dependence of the average distance 6 A from the magnetization threshold X C (U) for a Poissonian (squares) and 
Wigner-Dyson (diamonds) orbital distributions. A C (Z7) is extracted from Fig. [l^ as the value at which P(a > 0) = 0.01 and 
\c{U) + S\ corresponds to P(a > 0) = 0.5. The linear fits do not extrapolate to zero since P(a > 0) = 0.5 requires a finite SX 
at U/A = (see Figs. |l|). 
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FIG. 19. Magnetization parameter A vs. disorder strength W/V, for a Hubbard interaction and a two-dimensional 10x10 
(full circles), 20x20 (empty squares), 50x50 (full diamonds) and 80x80 (empty triangles) Anderson lattice. One clearly 
differentiates three regimes : (I) At small disorder, A increases due to a crossover from ballistic to diffusive behavior (see text). 
(II) At intermediate disorder, exchange and fluctuations compensate each other so that A is size-independent. (Ill) At large 
disorder one-body states are strongly localized over very few sites, which kills the off-diagonal fluctuations faster than the 
exchange and the latter dominates again. 
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FIG. 20. Magnetization parameter A vs. disorder strength W/V, for a Hubbard interaction and a three-dimensional 6x6x6 
(full circles), 8x8x8 (empty squares), 12x12x12 (full diamonds) and 15x15x15 (empty triangles) Anderson lattice. One 
clearly differentiates three regimes : (I) At small disorder, A increases due to a crossover from ballistic to diffusive behavior (see 
text). (II) At intermediate disorder, fluctuations seem to take over and A decreases with system size. (Ill) At large disorder 
one-body states are strongly localized over very few sites, which kills the off-diagonal fluctuations faster than the exchange and 
the latter dominates again. 
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FIG. 21. Effect of the interaction range on the magnetization parameter A for a two-dimensional 15 x 15 Anderson lattice and 
ratio Uo/Ui = 1 (full circles), 4 (empty squares), 9 (full diamonds) and oo (Hubbard interaction - empty triangles). Increasing 
the interaction range leads to a stronger increase of the fluctuations than of the exchange, resulting in a lowering of A. 
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FIG. 22. Schematic of the conductance peaks in a 2D lateral quantum dot as a function of an in-plane magnetic field gfiBB/A 
for m = 14, XU/A = 2 (left) and 0.4 (right) and U/A = 0.4 corresponding to the addition of the n = 3,4, ...9 th electrons 
(from bottom to top). The dashed lines indicate slopes of ±g/is/2 and serve as a guide to the eye. At weak exchange (right), 
electrons are piled up on the orbitals according to the Pauli prescription, so that the spin is always minimized and the peaks 
move in parallel to the dashed lines. The sequence of magnetization difference between consecutive peaks Sa z (n) — (— l) n /2 
results in peaks moving in opposite direction at low field. Larger slopes appear at larger exchange strength (left) and the 
successive addition of electrons of same spin results in conductance peaks evolving in parallel at low field. Note that due to 
the subtraction of the average charge-charge interaction, the model does not reproduce the charging energy so that the vertical 
distance between consecutive peaks is arbitrary. 
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FIG. 23. Representation of the a z = many-body yrast state with good total spin a at U/A — 0. The state consists of a 
filled fermi sea with n/2 — a doubly occupied orbital and a layer of 2a orbitals where particles are paired tripletwise. 
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FIG. 24. Transitions that do not exist for the cr^O yrast state. 
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